Hunting for CP violation with untagged charm decays 
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We construct a CP-violating observable which does not require flavor or CP tagging of the initial 
state. The proposed decay asymmetry could be measured at both threshold and non-threshold 
charm physics experiments and provide better sensitivity to small CP-violating parameters in charm 
decays. 



I. INTRODUCTION 

One of the most important motivations for stud- 
ies of CP violation in charm decays is the possi- 
bility of observing signals of new physics. It is 
expected that CP-violating amplitudes generated 
by Standard Model (SM) interactions are numer- 
ically small 0, so observation of CP-violation in 
the current round of experiments will provide a 
"smoking gun" signal of new physics, even if its 
source is not clearly identified. This robust predic- 
tion follows from the fact that weak decay of the 
charmed meson or baryon is governed by the real 
2x2 Cabbibo matrix, as all quarks which build 
up the hadronic states in non-lcptonic weak decay 
belong to the first two generations. The only pos- 
sible SM CP-violating amplitude comes from the 
virtual b quarks in penguin or box-type diagrams, 
but it is strongly suppressed by the the small com- 
bination of Cabbibo-Kobayashi-Maskawa (CKM) 
matrix elements V*ijVub ^ 0{\^), where A = 0.2 is 
the Wolfenstein parameter. This small "Standard 
Model background" makes charm decays a valu- 
able tool in searching for CP-violating effects of 
new physics, especially since the acquired datasets 
available in charm physics experiments are usually 
quite large. 

Let us first review the relevant formalism in or- 
der to collect necessary formulas. As in B-dccays, 
CP violating contributions can affect charm tran- 
sitions in three distinct ways 0. 

(1) CP violation can affect D decay amplitudes. 
This type of CP violation (also called "direct CP- 
violation" ) occurs when the absolute value of the 
decay amplitude Af for D to decay to a final state 
/ is different from the one of corresponding CP- 



conjugated amplitude, i.e. \Af \ 7^ \Af\. This type 
of CP-violation occurs in both charged and neu- 
tral D-decays and is induced by AC = 1 effective 
operators. The easiest way of observing it is by 
examining the rate asymmetries. 
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l + \Aj/Af\ 

where represents the D'^ — + i decay rate, while 
Ti represents the — > i rate. A two-component 
decay amplitude with weak and strong phase dif- 
ferences is required for this type of CP violation. 



Af =Af 



(2) 



Here Af corresponds to the strong phase differ- 
ence and 9 corresponds to the weak phase differ- 
ence between the CP-conserving (-4/) and CP- 
violating parts of the decay amplitude and ry 
represents the (small) ratio of their absolute val- 
ues. Note that Eq. (0) is a scale and scheme- 
independent way to write a non-leptonic decay 
amplitude. While no reliable model-independent 
predictions exist for the Af, it is believed that it 
could be quite large due to the abundance of light- 
quark resonances in the vicinity of the Z3-meson 
mass inducing large final-state interaction (FSI) 
phases. As the most optimistic model-dependent 
estimates put the Standard Model predictions for 
the asymmetry Acp < 0.1% 0, an observation 
of any CP-violating signal in the current round of 
experiments will be a sign of new physics. Cur- 
rent FOCUS, CLEO, and Belle/BaBar results put 
rather stringent bounds on Acp ■ For example, for 
a state K^K^ the direct CP asymmetry is U 

AcpiK+R-) ^ (0.0 ±2.2 ±0.8)%. (3) 
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Searches for this type of CP-violation in neutral 
charm decays must deal with the tagged data sam- 
ple, which means that only Z?-mesons tagged at 
production, typically in D* — > D^n decay, could 
be used in this analysis. This imposes restrictions 
on the available dataset of D's which could be 
used for this analysis. 

(2) CP violation can affect the D'^ — mixing 
matrix. This type of CP violation occurs when 
effective operators that change D'^ into D'^ states, 
i.e. generate the mass and width splittings for the 
mass eigenstates of the D° — £)° system, have CP- 
violating coefhcients. This results in the mixing 
of flavor eigenstates into the mass eigenstates, 



width differences of Di and -D2, 



\D.)=p\D°)±q\D') 



(4) 



Here the complex parameters p and q are the 
off-diagonal elements in the phenomenological 
parametrization of the — mass matrix ^ 
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CP violation in the mixing matrix occurs when 
2 
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It is convenient to define two dimcnsionless vari- 
ables X and y which are the normalized mass and 



m2 — mi 



2T 



(7) 



where mi, 2 and ri,2 are the masses and widths 
of the corresponding mass eigenstates, i.e. the 
eigenvalues of the mixing matrix in Eq. ISJ and 

r = (ri + r2)/2. 

This type of CP violation in charm can be ob- 
served most cleanly by looking for semileptonic 
decay asymmetries Asl, like the one in Eq.lQJ 
with / = Xtv. It is easy to see that Asl = 
{I - Rl,) I {I + Rl^) 0. This asymmetry is ex- 
pected to be tiny in both the SM and many of its 
extensions. 

(3) CP violation can occur in the interference 
of decays with and without mixing. This type 
of CP violation is possible for a subset of final 
states to which both 13° and can decay. It is 
usually associated with the relative phase between 
mixing and decay contributions. It can be studied 
by examining the time-dependent version of rate 
asymmetry Acp{f) of Eq. 

Studies of — oscillations offer a conve- 
nient probe of CP violation in the charm system. 
Using the notations of let us write the time- 
dependent decay rates of D° and to a given 
final state /. Since x,y <^ 1 0,13 and denoting 

r = rt, 



J 



Tf{t)^Rl\Af\\-^ 



Tf{t)^R;^'\Af\'e 



2 „-r 



{x' + y') - ix' - y') 



+ [Re (A/) y + Im (Xj) x]T + {x' + y') - {x' - y') |A/ 1 



^7 
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(8) 



where for a given final state /, CP violation is 
parametrized by 

Xj'^P^ = VRR-!e-^^'^^^+0{rf), 
qAf 

^7 = ^ VRR,ne-'('"^^ + 0{rj), (9) 

/ pAj f 

where Af and Af are the amplitudes for D'^ f 
and D'^ f transitions respectively, S is the 
strong phase difference between Af and Af, and 



i? = I ^/ / A/ 1 . Here (j) represents the convention- 
independent CP-violating weak phase difference 
between the ratio of decay amplitudes and the 
mixing matrix. The corresponding expressions 
for rj(i) and ^jit) can be found by substituting 
/^/in Eqs. ®. 

Eqs. © give the most general expression for 
time-dependent decay rate up to 0{x'^, y"^). How- 
ever, parameters of these equations (i.e. coeffi- 
cients of X and y to the appropriate power) scale 
differently for singly Cabbibo suppressed (SCS) 
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or doubly Cabbibo suppressed (DCS) decays. For 
instance, R - 0{X'^) for DCS decays like L»° ^ 
K~^Tr~, while R ^ 0{1) for SCS transitions such 
as TT+p'. This implies that in the stud- 

ies of a particular DCS or SCS transition, some 
of the terms in Eqs. ^ could be neglected. CP- 
violating parameters could be extracted by com- 
paring the time-dependent rates of D'^ and 
decays in Eq. ||SJ). 



violating parameters. Let's concentrate on the de- 
cays of D-mesons to final states that are common 
for D'^ and If the initial state is not tagged 
the quantities that one can easily measure are the 
sums 



r,(t) + r,(i) 



(10) 



for i = f and /. A CP-odd observable which can 
be formed out of is the asymmetry 



II. UNTAGGED SIGNALS OF CP 
VIOLATION 

The existing experimental constraints 0| 
demonstrate that CP-violating parameters are 
quite small in the charm sector, regardless of 
whether they are produced by the Standard Model 
mechanisms or by some new physics contributions. 
In that respect, it is important to maximally ex- 
ploit the available statistics. It is easy to see that 
the rate asymmetries of Eq. require tagging of 
the initial state with the consequent reduction of 
the dataset. 

Another way of looking for CP-violation in- 
volves methods which are being discussed in con- 
nection with CLEO-c tau-charm factory measure- 
ments. They do not require initial state flavor 
tagging but rely on the fact that at the threshold 
charm factory initial — D° state is prepared in 
the state with definite CP. An observation of a fi- 
nal state of the opposite CP would automatically 
imply CP-violation. These signals were discussed 
in . Since they involve CP- violating decay rates, 
these observables are of second order in the small 
CP-violating parameters, a challenging measure- 
ment. 

We propose a method that both does not re- 
quire flavor or CP-tagging of the initial state and 
results in the observable that is first order in CP 



CP 



f-^j _ N{t) 



+ Sj D{t) ■ 



(11) 



We shall consider both time-dependent and time- 
integrated versions of the asymmetry Hll|) . Note 
that this asymmetry does not require quantum co- 
herence of the initial state and therefore is accessi- 
ble in any D-physics experiment. From Eq. ((SJ it 
is expected that the numerator and denominator 
of Eq. (|lll) would have the form. 



N{t)=J:f - Ej = [A + BT + CT^] 



D{t)= e- 



\Af\ 



^7 



^7 



At 



.(12) 



Integrating the numerator and denominator of 
Eq. I|ll|l over time yields 



A^p{f)^-[A + B + 2C], (13) 

where D = r dt D{t). 

Both time-dependent and time-integrated 
asymmetries depend on the same parameters 
A, _B, and C . Since CP-violation in the mixing 
matrix is expected to be small, we follow |^ and 
expand R^^ = 1 ± Am. The result is 



A= [\Af\''- 
- 2yVR 



A- 



A-. 



1 - 



A-. 



B 







2\ / 2 
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J — cos(/)cosi5 1 \Af\ — 


^7 
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/\Afr]+Rl^l- Aj /\Af\ 

+ 0{AmX,rfX,...), (14) 



A- 



^7 



\AA' 



= —A + 0{Arnx'^,Arny^). 



We neglect small corrections of the order of 0{AmX,rfX, ...) and higher. It follows that 
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Eq. (|14|l receives contributions from both direct 
and indirect CP-violating amplitudes. Those con- 
tributions have different time dependence and can 
be separated either by time-dependent analysis of 
Eq. Hll|l or by the "designer" choice of the final 
state. Note that this asymmetry is manifestly first 
order in CP-violating parameters. 

In Eq. (|14|1 . non-zero value of the coefficient A 
is an indication of direct CP violation. This term 
might be important for SCS decays. The coefh- 
cient B gives a combination of a contribution of 
CP violation in the interference of the decays with 
and without mixing (first term) and direct CP vi- 
olation (second term). Those contributions can 
be separated by considering DCS decays, such as 
D K'-*K or D K''*'^p, where direct CP vio- 
lation is not expected to enter. The coefficient C 
represents a contribution of CP-violation in the 
decay amplitudes after mixing. It is negligibly 
small in the SM and all models of new physics 
constrained by the experimental data. Note that 
the effect of CP-violation in the mixing matrix on 
A, B, and C is always subleading. 

Eq. p4|l is completely general and is true for 
both DCS and SCS transitions. For an experi- 
mentally interesting DCS decay Z?" K~^Tr~ we 
can neglect direct CP violation and obtain a much 
simpler expression, 

A = 0, C = 

B= - 2y sin (5 sin (/) |l'x+,r-|^- (15) 

This asymmetry is zero in the flavor SU{3) f sym- 
metry limit, where 5 = j^. Since SU{i)F is 
badly broken in Z3-decays, large values of sin 5 0| 
are possible. At any rate, regardless of the theo- 
retical estimates, this strong phase could be mea- 
sured at CLEO-c ^12,]. It is also easy to obtain 
the time-integrated asymmetry for Ktt. Neglect- 
ing small subleading terms of C(A'*) in both nu- 
merator and denominator we obtain 

Agp(X7r) = -?/sin(5sin0/R. (16) 

It is important to note that both time-dependent 
and time-integrated asymmetries of Eqs. (|15|l and 
(|16|l are independent of predictions of hadronic pa- 
rameters, as both 5 and R are experimentally de- 
termined quantities and could be used for model- 
independent extraction of CP-viola ting phase (j). 
Assuming R ^ 0.4% and 5 40° [ll|| and y ^ 
1% one obtains \A^p{Ktt)\ - (0.04%) sin0. 
Thus, one possible challenge of the analysis of the 
asymmetry Eq. (|16|l . is that it involves a differ- 
ence of two large rates, E^+tt- ^-nd Y,k-tt+ i which 



should be measured with the sufficient precision 
to be sensitive to A^p, a problem tackled in de- 
terminations of tagged asymmetries in I? ^ Kn 
transitions. 

Alternatively, one can study SCS modes, where 
i? '--^ 1, so the resulting asymmetry could be 
0(1%) sin (/). However, the final states must be 
chosen such that A^p is not trivially zero. For 
example, decays of D into the final states that 
are CP-eigenstates would result in zero asymme- 
try (as Tf = Tj for those final states) while decays 

to final states like K^K*~ or p+7r~ would not. 

The final state / can also be a multiparticle 
state. In that case, more untagged CP-violating 
observables could be constructed. For instance, 
three body decays can exhibit CP-violating Dalitz 
plot asymmetries, like the i?+ ^ asym- 
metry of the Dalitz plot in the decay D 
KsTT^ {E+ )tt~ {E- ) . Similar studies of Dahtz plot 
asymmetries in B^-decays were suggested in [l^l ^ • 
Untagged studies of Dalitz plot population asym- 
metries resulting from the enantiometric interme- 
diate states were proposed in |lj| to study direct 
CP-violation in Bd decays. The use of untagged 
samples were also proposed to measure AF and 
CKM phase 7 in Bg decays ^, where large 
values of Xs cause these terms to oscillate rapidly 
and, effectively to cancel out for sufficiently large 
t. The situation in charm transitions is exactly 
opposite. 

As any rate asymmetry, Eq. Hll|l requires ei- 
ther a "symmetric" production of D'^ and , a 
condition which is automatically satisfied by all 
pp and e^e" colliders, or a correction for D'^ / 
production asymmetry. 



III. CONCLUSIONS 

We propose a method of searching for CP- 
violation in charm decays which does not require 
either fiavor or CP tagging of the initial state. The 
resulting asymmetry is first order in CP-violating 
parameters, which is important for charm tran- 
sitions. The unique feature of the asymmetry of 
Eq. (|11|) is that, apart from the small CP-violating 
phase, it could be sizable even for two body final 



^ That study, however, neglected the contribution due to 
Ufetime difference y (which a good approximation in Bj^ 
decays), leading for D-meson transitions. 

^ I thank Y. Grossman for pointing this reference to me. 
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states. This occurs because of the large SU{3)f 
symmetry breaking in charm transitions. 
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